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Abstract 

Wc show that the late time Hartle-Hawking wave function for a free massless scalar in a fixed 
de Sitter background encodes a sharp ultrametric structure for the standard Euclidean distance 
on the space of field configurations. This implies a hierarchical, tree-like organization of the state 
space, reflecting its genesis as a branched diffusion process. An equivalent mathematical structure 
organizes the state space of the Sherrington-Kirkpatrick model of a spin glass. 



1 X+sosceles Inception 

As envisioned in [1], an exponentially expanding universe, such as our own [2, 3, 4, 5], continuously 
breaks apart: its wave function constantly separates in new branches in which super-horizon modes 
of the metric and other light fields take on definite expectation values. Cluster decomposition (fac- 
torization of correlators at large spatial separation) holds within each separate branch but not for 
the full quantum state obtained by superimposing them. This process can be effectively thought of 
as a branched diffusion process [6] , in which quantum fluctuations of light fields in a given branch get 
stretched out to super-horizon scales, freeze and classicalize, while falling out of causal contact with 
each other. This spawns an offspring of new branches, which in turn give birth to new branches, and 
so on. 

The result — pictorially at least — is a tree-like structure on the state space of quantum fields 
on inflating spacetimes, and in particular on de Sitter space [7, 8, 9]. This should be distinguished 
from the tree-like causal structure of de Sitter space itself or models thereof [10, 11], where causally 
disconnected patches of space continuously spawn new causally disconnected patches of space. The 
two are of course related, as they both are produced by the peculiar generative dynamics of inflation. 

The dynamical tree structure should be encoded somehow in the wave function of the universe 
at late times, and thus, assuming some version of a dS/CFT correspondence [12, 13, 14, 15], in 
the partition function of the putative dual field theory at X^. In particular it should emerge from 
the Hartle-Hawking wave function of for instance a massless free scalar on de Sitter space. At first 
sight this seems hard to imagine, as the Bunch-Davies/Hartle-Hawking vacuum wave function of a 
massless free scalar is an almost trivial Gaussian at all times [16, 17]. Nevertheless, we will show in 
this paper that it emerges in a very sharp sense. 

The approach we take is inspired by methods developed to analyze the state space of spin glasses 
[18], in particular the remarkable Parisi solution [19, 20, 21] of the Sherrington-Kirkpatrick model 
[22]. A crucial element in this solution was the presence of ultrametricity [21, 23] in the state 
space: any three "pure" equilibrium states satisfy the property that their distances (measuring 
dissimilarity in local magnetization) form an isosceles triangle, with the unequal one the shortest 
of the three. Equivalently, these states can be organized as the leaves of a tree, with the distances 
corresponding to how far back on the tree one has to go to find a common branch. Yet another 
way of phrasing this is that the state space forms a hierarchy of nested clusters, with distances 
given by the size of the smallest common cluster. The Parisi solution involves an auxiliary flow 
equation reminiscent of a Hamilton-Jacobi equation in an expanding space or an RG equation (as 
reviewed in [7]), which generates the ultrametric structure. Although it is computationally intractable 
to enumerate the equilibrium states explicitly, it is possible to compute exact n-point probability 
distributions of distances between equilibrium states, sampled from the Boltzmann-Gibbs measure. 
Such distance distributions serve as order parameters for the spin glass phase. Ultrametricity shows 
up in the fact that the three point distance distribution has support entirely on isosceles triangles. 
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We will compute the one point and three point Euclidean field space distance distributions for 
a free massless scalar in dS^+i for arbitrary d, directly from the Hartle-Hawking wave function 
evaluated at late times, and without making any assumptions about effective classical behavior. In 
the Minkowski case, after fixing the zeromode, the distance distribution becomes a delta-function 
peaked at zero, since there is a unique vacuum, which satisfies cluster decomposition. But in de 
Sitter, as we will see, these distributions acquire a finite width, with a scale set by the de Sitter 
temperature. This confirms the Bunch-Davies vacuum does not satisfy cluster decomposition, but 
that it decomposes instead in an infinite set of branches which do. (See e.g. [7] for a more detailed 
explanation of the relation to the clustering property.) The distributions are always of Gumbel type. 
Gumbel distributions arise as extreme value distributions, i.e. distributions of maxima or minima 
of a set of random variables. They pop up in a wild variety of physical systems (see for example 
[24]). Interestingly, the same Gumbel distribution was found very recently in [25] as the distance 
distribution for a free massless scalar in dS2, but for an a priori rather different distance measure, 
suggesting it is a universal feature. 

More strikingly, we will demonstrate directly that the three point distance distribution exhibits 
nontrivial ultrametricity: Given two distances di2 and d23 with di2 less than d23 by more than a few 
dS units, the conditional probability for the third distance d^i peaks sharply on ^23, the longer one 
of the two. The dynamical origin of this hidden tree structure in the wave function at a fixed time 
is simple: The fields drift away from each other at a steady rate, so the distance between different 
states (labeled by large scale field vevs) at some fixed time slice is proportional to how far back in 
time one has to go to get to a common ancestor. What is remarkable is that this can be seen in a 
simple Gaussian wave functional. 

2 Defining Distance 

We begin with a brief discussion of the notion of distance between field configurations. We will 
restrict to spatially flat FLRW cosmologies. 

2.1 Geometry 

The (d + l)-dimensional FLRW geometry with scale factor a{r]) is given by: 

ds"^ = a{r]f {-dri^ + jijdx'dx^) , i = l,2,...,d . (2.1) 

Constant conformal time rj slices are given by three-dimensional slices of constant curvature with 
induced metric a(?7)^7jj. We will restrict our analysis to jij = 5ij in what follows. Furthermore we 
assume periodic boundary conditions Xi ^ Xi + L. The pure de Sitter geometry with flat slicing, 
which is our main focus, has a(ry) = £/rj with ij G [— oo,0]. Future inflnity X"*" resides at r/ = 0. The 
quantity £ is the de Sitter length which is related to the cosmological constant A as £^ ~ 1/A. 
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We consider the evolution of non-interacting massless scalar fields in this background geometry. 
Massless fields experience quantum fluctuations which freeze out and reach all the way into the 
future. This leads to a distribution of possible field configurations on the late time spacelike slice. 
The massless scalar is the simplest case to consider, and it adequately models the behavior of metric 
perturbations. 

2.2 Distances in field space 

Given two field configurations ipi (x) and ip2 (x) on a constant time slice, we define a distance between 
them as follows. First we define new fields (pi and 02 by subtracting the zeromode and coarse graining 
over some finite scale by convolution with some window function. The latter is equivalent to imposing 
a physical size UV cutoff on the field modes. Thus: 

0ix) = [ dS w{y) ^{x + y)-^ [ d'^x ip{x) , (2.2) 



where w{y) is a window function (e.g. w{y) oc ' ) which we imagine to have a width of the order 
of the Hubble radius. The zeromode has no natural normalizable ground state and it is unobservable; 
subtracting it here will allow us to conveniently discard this mode altogether for our purposes. 
We then define a distance between two ^pi{x) and f2{x) as 

di2 = d[ipi, ^2] = ^ [ d'^x ((^i(f) - 02ix)f . (2.3) 



As it stands, in de Sitter space, the expectation value of di2 diverges linearly with proper time 
in the far future. To avoid this issue, we further subtract off the average (c?i2) with respect to the 
probablility distribution dictated by the quantum state. Thus we define the following "renormalized" 
distance: 

S12 = di2 - {di2) . (2.4) 

Relatively similar configurations will have negative values of 5i2 whereas relatively dissimilar vacua 
will have positive values. 

3 Hartle-Hawking 

In this section we review the construction of the Hartle-Hawking vacuum wave function [17], which 
we take to be state our massless scalar field is in. 

3.1 Scalar Solutions 

A minimally coupled massless free scalar ip{r], x) = e^^'^ ipj^{r]) in an FRW background of the form 
(2.1) has action 
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where the prime denotes rj differentiation. The zero mode A; = drops out of the distance (2.4), which 
allows us to discard it altogether. In what follows it is always understood that /c = is excluded 
from sums and products. The general solution to the equations of motion for fj^iv) can be written 
as a linear combination of mode functions: 

^^{r,) = Apkiv) + AZvUv) ■ (3.2) 

In Minkowski space M^+i (a(??) = 1), the canonical mode functions are Vk (ry) = ^7=1 e^^^ while in 
dSd+i (a(r?) = -i/rj) they are 



vUn) = fr^'^-'^^' {-If' H^,%ik,) , (3.3) 

where H^^\y) is the Hankel function of the first kind. Asymptotically Hn\y) ~ when 
y — — oo. Specifically in dS2 this becomes Vk{ri) = ^^^^ e^^^, and in dS4 

The normalization is the natural one for canonical quantization but will actually not matter for our 
purposes. 

3.2 Vacuum state 

The standard Bunch-Davies vacuum state [16] of a free scalar in dS^+i, evolved to a time rj = rjQ, 
can be represented as a Gaussian Hartle-Hawking wave functional: 

^-ir/oi^] oce^^t*'"^! , (3.5) 

where is a solution to the Hamilton-Jacobi equations, obtained from the action (3.1) evaluated 
on the complex solutions 

<^fc(^) = -^4-^ ' V<Vo- (3.6) 
Vk{rio) 

This satisfies the boundary conditions <^g(??o) = </>^ and f^iv) ~ e*'^'' for rj — )• — oo. Integrating (3.1) 
by parts reduces the on-shell action to a boundary term and identifies 

W[wA] = ^a(r?o)'^-^L'^X;(l°g^'^)'(^o)l'^A;l' ■ (3-7) 

it 

Indeed for this choice the wave functional e*^ is a properly damped Gaussian, mimicking for each 
mode k the canonical Minkwoski vacuum in the far past. 

The probability functional of a given configuration (p at some time r/ = 770 is thus given by 

V[<P] = m<P]\^ oc e-2EE&l0gP ^ ^ Re ['-a''-^ L^^logvk)'] . (3.8) 
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For example in dS2 we have /3fc = kL/2 and in dS4 we get /3a,. = ^^(LA;)^/2(1 + A;^?7q). More generally 
in the late time limit r/o — we find for dS^+i 

h = \{Lk)H''-\ (3.9) 

4 Distance Distributions 

We are now ready to compute the distance distributions. We begin by making more precise the 
notion of "branches" of the wave function of the universe, as states characterized by a definite value 
for the coarse grained fields. We point out the analogy with thermodynamic pure states in statistical 
mechanics, in particular spin glass theory. Following this analogy, we define a notion of distance on 
the state space and introduce their probability distributions. In the theory of spin glasses these serve 
as order parameters capturing the spin glass phase. We proceed by computing these distributions, 
first for dS2, for which exact results can be derived, and then for the general dS^+i case. The 
resulting three point function exhibits ultrametricity. 

4.1 State space and replica reasoning 

Let be the coarse grained field of (2.2). In the Bunch-Davies/Hartle-Hawking vacuum, equal time 
correlators {(p{x)ip{y)) are logarithmic and do not vanish at large separation — in other words the 
Bunch-Davies vacuum does not satisfy cluster decomposition. Related to this is the fact that (^(x) 
does not have a definite value: starting from the Bunch-Davies vacuum, it could freeze into a huge 
number of possible profiles at late times. 

This is analogous to how the usual Boltzmann-Gibbs state does not satisfy cluster decomposition 
and order parameters do not have definite expectation values whenever there exist different equilib- 
rium states (like the spin up/down states of the Ising model at low temperatures). In analogy to 
how in spin glass theory one decomposes the Boltzmann-Gibbs state Vbg ~ e~^^ into "pure" states 
which do satisfy clustering [18], we imagine decomposing the Bunch-Davies state into such states. 
More precisely, we imagine decomposing the Hartle-Hawking probability functional V[(l)] (3.8) as 

V[(p] = ^WaVa[(l>], ^u;, = 1. (4.1) 

a a 

Here a is an abstract label for the states exhibiting the clustering property and definite expectation 
values {(p)a for the field ip. Happily, as in the spin glass case, it is in fact possible to extract detailed 
information about the structure of the state space without having to perform this decomposition 
explicitly. The key idea is to consider probability distributions of distances or overlaps between pairs 
of pure states. In the case at hand, the distance between two states a and /3 is 

dafS = d[{ip)a,{y^)l3], (4.2) 



6 



where d[-, ■] is as defined in (2.3) and {tp{x))a is the local vev of ip{x) (or more precisely of (^(x)) in 
the state labeled by a. The probability distribution for finding a distance D between two states is 
then 

P{D) = Y,WaWp5{D-d^p), (4.3) 

a/3 

and similarly P{Di,D2, D^) = X^^^^ WaWpw^ S{Di — dap) 5{D2 — dp^) 8{D^ — d^a)- As they stand, 
these formal expressions cannot be evaluated. However, as in the spin glass case, in the thermody- 
namic limit, one can rewrite this as the distance distributions between configurations of a pair of 
replicas sampled directly from the Bunch-Davies vacuum state: 

P{D) = {5{D-di2))^ I P</.«P</.(2)|vI/HH(0('))Pl^HH(0('^)|'5(I)-d[</.W,0(2)]), (4.4) 



and similarly P{Di, D2, D^) = {5{Di — d23) 5{D2 — d3i)5{D3 — ^12)), now involving three replicas. 
As we will see, these distributions are easily computed. 

As mentioned already in section 2.2, the late time expectation value of di2 diverges, so we consider 
the renormalized distances 612 instead, as defined in (2.4), with corresponding probability distribution 
P(A) = {6{A-6i2)). 

4.2 Generating function 

We will focus on the distance distribution in the late time limit t/q — of de Sitter space. In the late 
time limit, the number of Hubble and coarse graining volumes goes to infinity, so this can also be 
thought of as the thermodynamic limit. In this limit, the UV cutoff on k imposed by coarse graining 
ip is pushed to infinity, so as long as we are computing UV finite quantities, we can remove the cutoff 
completely, such that only the zeromode subtraction remains in (2.2). 

In order to compute -P(A) it is convenient to first compute the moment generating function 
G{s) = (e^'"^^^), and from this obtain P(A) = ^ J1^^g^^G{s). For the probability functional 
(3.8), the function G{s) is given by a simple Gaussian functional integral: 

k 

The primed product runs over unordered pairs {k, —k) with k ^ 0, which factorizes the integration 
over the independent complex variables (pj^ = (jf ^. The resulting Gaussian integrals are proportional 

to _^ 

Z{k,s)^det(^^'^' ] (4.6) 

and the proper normalization is obtained by dividing by (1), i.e. A/fc = 1/Z{k,0). Thus, 
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Figure 4.1: Plot of the Gumbel distribution in equation (4.10). 



with the factor e^^^^ coming from e^**"'^^^. In a similar fashion, one obtains the moment generating 
function for the triple distance distribution: 

p(si+S2+S3)//3fc 



^g-S3<5l2g-S2<5l3g-Sl(523^ 



+ {Sl + S2 + S^)/(3k + 3(S2S3 + SlS2 + SiS3)/4/3| 



(4.8) 



4.3 dSa 

In the case of dS2 we have = kL/2 = vrn with n G Z"*", and the infinite product (4.5) can be 
expressed as follows: 



1 + 



TT 



(4.9) 



-'j^ 1 + s/vrn 

where = 0.577216... is Euler's 7 constant. The inverse Laplace transform of the above expression 
provides -P(A). This is easily performed using Euler's integral expression for the F function given 
by T[z\ = t^~^e~^dt. The resulting distribution is a Gumbel distribution: 



P(A) 



ds 
2Td 



g«A eW^r[l + s/7r] = 7re 



-(7rA+7E)-e-('^^+TB) 



(4.10) 



The Gumbel distribution is an extreme value distribution. It appears in the analysis of wide variety 
of physical systems as discussed e.g. in [24]. We depict it in figure 4.1. There is a clear asymmetry: 
It falls of much faster for negative values than positive ones. This reflects the fact that there are 
exponentially more dissimilar configurations then there are similar ones. 

The three point distance distribution can be similarly analyzed in dS2. Once again, the infinite 
product can be explicitly done. We find: 



G{s) = (e 



3<'12g-S2<5l3 



g-si523|) = g7i5i(s")Ap[]^ _^ L(s)+Q(s) jp^j^ _^ L{s)-Q{s) - 



2tt 



(4.11) 
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Figure 4.2: Left: Conditional probability P(A3|Ai,A2) for Ai = —2. Larger values of P corre- 
spond to lighter colors. The sharp peaking on A3 = max{Ai, A2} indicates ultrametricity. Right: 
P(A3|Ai, A2) for Ai = —2, —3, —3.5 from top panel to bottom panel and A2 = —2, —1,0, 1, 2 from 
left peak to right peak. The peaks get sharper for lower Aj. 

where we defined the permutation invariants: 

L{x) = Xi+X2 + X3, Q{x) = ^\{xi- X2Y + 5(2:2 - 2:3)2 + 5(X3 - xi)2 . (4.12) 

The three point distribution i-'(A) is the inverse Laplace transform of (4.11). We were unable to do 
this exactly. After replacing the L-functions by their integral representation, a straightforward but 
somewhat tedious saddle point evaluation yields 

P(A) ocexp(-|L(A')-2e-5^(^')cosh[|g(A')]) , (4.13) 



where A' = vrA + 7^ and 7^ = (7^;, 7s, 7s)- 

The conditional probability to find the configurations 0^^^ and (/)(2) ^t a renormalized distance A3 
given that <\Pi and 4>^^'^ are at a renormalized distance Ai and that (^^^^ and (/>(^^ are at a renormalized 
distance A2 is 
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We depict this distribution in fig. 4.2. A striking structure emerges: the distribution tends to peak 
very sharply on isosceles triangles — in other words, the state space is essentially ultrametric! 

This can be seen explicitly from (4.13). First note that (5(A') > ^^^|Ai — A2I, so that if |Ai — A2I 
gets larger than about one, the cosh function is well approximated by an exponential. The second 
term in the exponential (4.13) is then minimized by minimizing — L(A') +2(5(A'). The minimum of 
this term is e"'^™'^^^!'^^} g^j^^j jg reached at 

A3 = max{Ai,A2}, (4.15) 

i.e. on an isosceles triangle with the unequal side the shortest of the three. When min{Ai, A2} < 0, 
i.e. at least two configurations are closer than average, there will be superexponential suppression of 
deviations from (4.15), and the distribution becomes strongly ultrametric. This is also clearly visible 
in the figures. Note that in this regime the total probability is always exponentially small as well. 
This is because we are conditioning on the intrinsically rare event of sampling two relatively similar 
configurations. The ultrametric structure is for this reason easy to miss if one simply plots the total 
joint distribution including the region where it gets maximal. 



4.4 dSd+i 

For higher dimensional de Sitter, the product (4.7) cannot be computed in closed form. We can 
however approximate its logarithm by an integral, taking into account the momentum quantization 

log{.-*)»l/A^ (_^_,„,(l + _i_)) , (4.16) 

where for dSd+i in the late time limit rjQ — t- 0, according to (3.9), we have (3{k) = ^{Lk)'^£'^^^ . On 
the other hand the volume of a shell in momentum space is ujd d{k'^) where cod = {d/2)\ volume 
of the d-dimensional unit ball. Happily this implies the integral takes a universal form for all d: 



The chosen value of the IR cutoff /3o corresponds to /cq = 3- • This integral is easily computed, and 
we can evaluate the inverse Laplace transform again in saddle point approximation, yielding again a 
Gumbel distribution: 

P(A) ^ exp[-K(A' + e-^')] , (4.18) 

where 

''^Y = ^(^' A'^7 + ^(2^^)'^-^A, 7 = M^^) - V'(^) , (4.19) 

with ip the digamma function. The shift by 7 and normalization factors do not follow directly from 
the saddle point computation but are inferred from the requirements (1) = 1, (A) = 0. In the case 
d = 1, we have k = 1 and A' = 7^; + vrA so this reproduces the exact result (4.10). For d = 3, we 
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have K = 2'K l?> and 7 ~ 0.257336. Notice that the characteristic scale for the distance distribution is 
set precisely by the dS temperature T = l/27r£ [26]. 

Finally the conditional probability for A3 given Ai and A2 can be obtained by similar approxi- 
mations. We find: 



P(A) ocexp [-K(|L(A') + 2e-5^(^')cosh[|Q(A'; 



(4.20) 

with K and A' as in (4.19). This is of the same form as (4.13). We conclude that up to smoothing 
effects at a scale set by the de Sitter temperature, the state space of a massless scalar in de Sitter 
space is universally ultrametric. 
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